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Abstract : People arc familiar with quantum mechanical reflection and transmission coefli- 
cient. In all those cases corresponding potentials are usually assumed as of constant height and 
depth. For the cases of varying potential, corresponding reflection and transmission coefficients can 
be found out using WKB approximation method. But due to change of barrier height, reflection 
and transmission coefficients should be changed from point to point. Here we show the analytical 
expressions of the instantaneous reflection and transmission coefficients. Here as if we apply the 
WKB approximation at each point, so we call it as Instantaneous WKB method or IWKB method. 
Once we know the forward and backward wave amplitudes we can find out corresponding wave 
function by calculating Eiconal. For the case of analytically complicated potential corresponding 
differential equation seems to be unsolvable analytically. If the potentials are well behaved then 
obviously these could be replaced by functions of simple expression of same behaviour which can be 
integrated analytically and the equation is now possible to solve analytically. 
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I. INTRODUCTION 

In the cases of quantum mechanical barrier problem reflection and transmission coefficients can be calculated. In 
those cases with the interaction of the potential field, one fraction of incident particle (or wave) is reflected back 
to the infinity (reflection coefHcient) and other fraction transmits into other side or may tunnel through the barrier 
(transmission coefRcient). If the barrier is of constant height and depth, the reflection and transmission coefficient 
can be c^asily calculated [1]. In the case of varying potential, from the far away of potential field, rcficc;tion and 
transmission coefficient can be calculated using WKB approximation [1,2]. In this present paper we will show how the 
reflection and transmission of the particle changes in the cases of varying potential point to point. Due to variation of 
potential, in each point particle should feel different potential field so the fraction of particle (or wave) which reflects 
and transmits will be changing at each point. Here we will calculate the analytical space dependent expression of 
reflection and transmission coefficients which indicates transmission and reflection of the particle at any arbitrary 
point with respect to the transmission of the particle at the immediate previous point. So the local values of incident 
and reflected wave amplitude can be found. With this, we shall present the solution of the differential equation (i.e. 
Schrodinger equation), containing that space dependent reflected and transmitted amplitude. 

Secondly, in our daily problems of physics and mathematics there are several differential equations which can not 
be solved analytically. If we want to solve those equations we need numerical methods. If, however, the nature of the 
coefficients of dcirivative terms or inhomogencous term namely factor, due to presence of those factors corresponding 
differential equation to be unsolvable, are well behaved then using the methods which will be explained here the 
differential equation can be solved analytically using WKB approximation method. Actual reason of non-analyticity 
is that corresponding factors can not be integrated analytically. 

We can consider, a particle is moving in a potential field whose analytical form is complicated. Obviously corre- 
sponding Schrodinger equation will be non-analytical. Although the reflection and transmission coefficients for that 
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particle can be calculated but the wave function can not be analytically calculated exactly, because for this job we 
need to integrate analytically the wave vector k, which includes the complicated potential. 

In next section we will discuss how the space dependent reflection and transmission coefficient can be calculated for 
the case of varying potential. In §3, we will establish approximate analytical solution of a Schrodinger equation with 
the presence of potential with complicated analytical form. In §4 we will illustrate an example where our method can 
be applied for finding solution and reflection-transmission coefficients. In §5 we will make our conclusions. 

II. METHOD TO FIND LOCAL REFLECTION AND TRANSMISSION COEFFICIENTS 

We consider an one dimensional Schrodinger like differential equation as 

^+kHx)y = (1) 

where, k^{x) = E — V{x), E = total energy, V{x) = potential energy. (2) 
Here, the independent variable x is varying from — oo to oo, i.e., nothing but cartesian co-ordinate. 
Using WKB method the solution of this equation is [1,2] 

A B 

y{x) = --j=exp{iu) + -j=exp{-iu) (3) 

where, 

{x)dx, (4) 



A and B are constants of integration. 

The reflection and transmission coefficients calculated by WKB method means the coefficients are calculated at a 
certain point, that means from a certain point what fraction of particle is transmitted inside and what other fraction 
is reflected outside. If we calculate this transmission and reflection of the particle from different points, results should 
be different. So the probability of transmission and reflection of the particle from point to point will be different 
for the case of varying potential. So the constants A and B should be changed from point to point due to change 
of effective potential which is felt by the particle. Here we will give an analytical space dependent expression for 
transmission and reflection coefficient by which one can find, how the reflection and transmission probability change 
from point to point for the case of varying potential. 

It is assumed that far away from the potential field (at x = oo), the potential barrier height is almost constant, so 
the refiection and transmission coefficients are almost same. In those regions we safely can choose A and B as pure 
constants. Again, all along the sum of transmission and reflection coefficients should be unity, so far away from the 
potential field the relations between A and B are, 

A = B + c{c = constant), (5) 
A2 + £2 = (6) 

It is very easy to check from incident and reflected current that transmission and reflection coefficients are ^ and 

respectively. 
Solving these two equation we get. 



c y/2k - c2 



A{x) = - + ^-^ (7a) 



B{x) = -- + (76) 

These A{x) and B{x) are very slowly varying functions at large x. These are indicating, at far away from the 
potential barrier (where potential is varying in very slow manner), how constants A and B are slowly changing. So 
from these slowly varying A and B, one can find out how transmission and reflection coefficients are varying at far 
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distance. Now, constant c can be calculated from the boundary conditions. From pure WKB approximation constant 
c can be calculated. From WKB method, far away from the potential field reflection and transmission coefficients 
and corresponding reflected and transmitted amplitude can be calculated. So constant c is immediately determined. 
Obviously, at a; = — oo these expressions of A and B will not be valid because of simplified assumptions. At the 
minimum value of x, i.e., x = — oo, we have to put another condition. 

It is clear that the potential field, which is felt by the particle, extended in between — oo to oo and assumed that 
at a; = — oo barrier height goes down to zero or constant value i.e., the potential is asymptotically flat. So the 
particle which reaches at x = — oo, should feel free. As a result, there transmission should be 100% and corresponding 
reflection is zero (hereafter, this will be called as inner boundary condition). By introducing this inner boundary 
condition in eqn. (7) we get, 

C{x) = ci + A{x) = ci + I + , (8a) 



r \/9k — 

D{x) = C2 + B{x) = C2 - 2 + 2 ^^^^ 

Here, constants Ci and C2 are introduced to modify the reflection and transmission coefflcients according to the inner 
boundary condition. 

Here, one necessary condition is the sum of reflection and transmission coefficients should be one. By modifying 
the coefficients, it is seen that 



C^{x) + D\x) = (ci + I) + (c2 - I) + (ci + C2)V2k^ + ^-^Y^ = h{x)[say] 
So, it is advisable to choose the modified coefficients of the wave function as follows: 



h{x) 



k{x) 



k{x). (9) 



a{x) = (10a) 



b{x) = (Wb) 



so that 

a'^{x)+b'^{x)=k{x) (11) 
and the transmission and reflection coefficients are and respectively which are explicitly written as: 

Now we describe, how to determine the constants Ci and C2. First of all to satisfy B{x = — oo) = (inner boundary 
condition) C2 is determined and fixed. Then introducing this C2 in h{x) and equating the expression of modified 
reflection coefficient [b'^ (x) / k{x)] and reflection coefficient valid only at x = oo [B'^{x)/k{x), which is correct at far 
away from the black hole] we get the other constant ci . These expressions for transmission and reflection coefl&cients 
are valid in whole region. The space dependent coefficients of incident and reflected waves, a{x) and b{x) are valid 
for arbitrary x. 

So the flnal form of the solution using WKB method is 

y(x) = — j^exp(iu) H j^exp(—tu). (13) 

Vfc Vk 

Here one important thing is to be noted that, WKB approximation is valid only if << fc- So one can calculate 
the space dependent reflection and transmission coefficient only if potential vary in such a manner that k satisfies the 
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above condition. Otherwise our method to find out space dependent amphtudes and coefficients is not valid because 
our method is based on WKB approximation. Equation (13) is valid at any point. It is indicating the solution and 
corresponding reflection and transmission coefficients for an arbitrary x as if at each point WKB approximation method 
is applied instantaneously. So our modified WKB method can be called as Instantaneous WKB Method or IWKB 
method [3]. For the cases, where << ^ is not satisfied, the potential can be replaced by a large number of square 
steps and space dependent refiection and transmission coefficients can be found [3]. Here, corresponding problem is 
reduced to as simple quantum mechanical barrier problem with multiple steps where at each step junction the wave 
functions and its derivatives of two separate regions should be continuous. Detail discussion is in Mukhopadhyay & 
Chakrabarti [3]. 



III. SOLUTION OF THE EQUATION 

In the previous section, although we have given the form of the solution with space dependent incident and reflected 

coefficients but the Eiconal u still yet to be determined. If the analytical form of the potential and corresponding 
wave vector k is well integrable analytically then u can be calculated immediately but sometimes the form of k may 
be complicated such that it can not be integrated analytically then u can not be evaluated unless very special cases 
are chosen. In this section we want to give an approximate solution of eqn. (1) using WKB (actually IWKB, as 
mentioned in last section) approximation even if the analytical form of k is not well integrable. As we mentioned in 
INTRODUCTION that if k'^{x) is well behaved then we can give an analytical form of the eqn. (1) using the following 
trick. If the nature of the coefficients of derivative and zeroth order derivative terms are well behaved but analytically 
complicated then those analytical form of the coefficients can be replaced by piecewise continuous analytical function 
with simple form and of same behaviour as exact complicated functions. After this replacement k will be of simple 
form and we can integrate it to find out u. Once we know the analytical form of u, we can present the analytical 
solution of the equation. 

In different ways we can replace the complicated functions. Any function can be expanded in terms of polynomial. 
It is seen that if the corresponding polynomial is of order two then the corresponding k can be integrated easily. But 
using one polynomial of order two, in the total range, the replacement of the complicated function can not be done. 
So in different regions, different coefficients of the polynomial are chosen i.e., different second order polynomials are 
chosen such that at the boundaries (where the polynomials change), value of the two different polynomials and its 
derivatives are same. In this way, by using many number of piecewise polynomials the analytically complicated well 
behaved functions can be replaced. 

If in the differential eqn. (1), the analytical form of V and corresponding k is of non-analytical form (complicated 
analytical expression) then potential V can be replaced by piecewise continuous polynomial as 



Vi{x) = ai+ bix + cix^ 



(14) 



and the corresponding k{x) is defined as 



kf{x) ^ {E - ai) - hix - cix^ 



(15) 



Integrating this we can find out analytical form of u as 



ui{x) = -2ciki - ci\/{E - ai)log 



a, 



ki + ^{E - ai) 



constant 



(16) 



Where index / in the coefficients of the polynomial is indicating that in different ranges of x different values of a, 6, 
c may be chosen. This type of replacement also can be done by some other well integrable functions according to 
the nature of the complicated functions. In the next section we will give an example, where, to solve the differential 
equation with the potential of complicated analytical form, we replace the complicated function (potential) by the 
function of simple analytical form otherthan second order polynomials. 

In this way we can give the analytical form of the solution of Schrodinger equation like equation (3) even if the 
analytical form of the potential is complicated. Since we have used the WKB approximation method (more clearly 
IWKB method), our solution only will be valid for ^ ^ << fc. If the potential vary in such a way that above condition 
is not being satisfied then the solution of the Schrodinger or Schrodinger like equation is not possible by this method. 
For those cases we have to employ some other methods, one of them is indicated in last section (by reducing the 
problem to simple quantum mechanical barrier problem). Sometimes for those cases, the Schrodinger equation can 
be reduced into Bessel equation and solution is possible by expressing in terms of Airy functions [1] . 
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IV. ILLUSTRATIVE EXAMPLE 



Dirac equation in Kerr geometry can be separated into radial and angular parts [4-5] . Then corresponding decoupled 
radial equation on some transformation of independent variable can be reduced into Schrodinger like equation [5] as 

+ (a^ - V{.)) = 0. (17) 

In comparison with Schrodinger equation with H = c = G = 1, wc can say V{x) is nothing but potential and is 
proportional to energy of the particle. Here the analytical form of the potential is complicated as 

where a — Kerr parameter, rrip — mass of the incident particle, a = frequency of the incident particle, M — mass of 
the black hole, 'X = separation constant, r = radial co-ordinate, A = — 2Mr + a^, w = r'^ + + am/a, x = f{r). 
For details see [5]. 

Now, although the analytical form of the potential is complicated but for particular set of physical parameter (such 
as, a, a, nip etc.) the nature of V is well behaved. So wc can replace it by suitable peace-wise continuous analytical 
function of simple form. We can choose it as second order polynomial as explained in previous section. But for 
convenience we choose here in different form. 

To show the nature of the potential, wc choose as an example one set of physical parameter, such as: 
a = 0.5; M = 1; rup = 0.8; orbital quantum number, I = i; azimuthal quantum number, m = — i; a = 0.8; A = 0.92 
from [6]. 

The simple analytical form which is chosen to replace the complicated form of the potential for this set of physical 

parameter is given as 



ci 



Vi{x)=ai + biexp[--] (19) 



and corresponding k and eiconal u are given as 



ki{x) = J{a^ - ai) - biexpi--), (20) 

ci 



u{x) = -2ciki{x) +ciy/{a'^ - ai)log 



VW~<n) + kix) 



a/(o-2 - ai) - ki{x) 



■ constant. 



(21) 



The behaviour of the functions of eqn. (18) and (19) are same. In case of (19) in different ranges of x different values 

of a, 6, c arc chosen. 

For these physical parameters, coefficients a, 6, c of mapping function in different ranges of x are given as: 

a = 0,b = -0.187354, c = -3.75 for x -co to 0, 

a = 0.603, b = 0.415646, c = 8.79 for x to 30, 

a = 0.629, b = 0.12690038, c = 26.3 for x 30 to 109, 

a = 0.63543098, /; = 0.037193439, c = 73.5 for x 109 to 208, 

a = 0.63543098, b = 0.2228925, c = 45 for x 208 to 310, 

etc. 

So finally wc can say the analytical form of the solution of the Schrodinger equation is possible for this set of physical 
parameter. Similarly for other sets of physical parameter one can find out solution following same method (by mapping 
the complicated analytical form of the potential by the simple form which may be polynomial or analytical function 
like eqn. (19) and which are analytically integrablc). 

The space dependent transmission and reflection coefhcients can be calculated following the method explained 
earlier. Now we will calculate the constants c, Ci , C2 by imposing boundary conditions for the potential of given set of 
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physical parameter. Using general WKB approximation method, from the far away of the black hole (say at x = 310) 
the reflection and transmission coefficients of the particle for the potential shown in Fig. 1 can be calculated as 
T = 0.299 and R = 0.701. Using these values of the coefficient constant c can be calculated as 

c = \/Tfc - \/M = -0.0913 (22) 

where k is the value of wave number aX x = 310. 

Using this c when we calculate the transmitted and reflected amplitude from (7a-b) (which are valid at large 
distance from the black hole) obviously it will not satisfy the inner boundary condition. To satisfy the inner boundary 
condition first C2 is introduced to reduce B{x) to zero at a; = — oo as C2 = —0.6765. Then as explained in section II, 
by matching the reflection amplitudes of eqn. (7b) and (10b) at infinity i.e., 

b{x) = B{x), (23) 

Ci can be evaluated as ci = 0.1639. 

Now all the constants are known. In Fig. 2 variation of instantaneous reflection and transmission coefficients 

are shown. It is seen that as the potential barrier lower down the transmission probability increases and reflection 
probability decreases. For the cases of other set of physical parameters where the potential barrier may be of different 
type, using the same method constants can be calculated. In the final solution with the insertion of the value of 
constants it can be checked that WKB approximation is still valid except in the regions where is close to V and 
this modified WKB (IWKB) solution satisfies the original differential equation if << fc. 

V. CONCLUSIONS 

Here, we show how reflection and transmission coefficients vary for the particle moving in spatially varying potential 
field. Actually we solve the Quantum Mechanical barrier problem, where the barrier height changes point to point. 
For solution, we introduce modified WKB approximation method namely instantaneous WKB method i.e., IWKB 
method where at each point, instantaneously, the WKB approximation method can be applied and corresponding 
reflection and transmission coefficients can be calculated. Since the potential changes point to point corresponding 
reflection and transmission coefficients change. We also indicate how Schrodinger equation is solved analytically where 
coefficients of derivative terms are of complicated analytical form. It is shown that if coefficients are well behaved 
then the the analytical solution can be set up by suitably mapping the coefficient-functions to fiuictions of simple 
form. In this way, any such kind of differential equation can have analytical form of the solution which could be useful 
for further study. 
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FIGURE CAPTION 

Fig. 1 : Variation of potential V (solid curve) with respect to radial distance x and total energy E (dashed curve) of particle. 
Fig. 2 : Variation of transmission (solid curve) and reflection coefficients (dashed curve) with respect to radial distance x. 
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